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Abstract. - We assume that markovian dynamics on a finite graph enjoys a gauge symmetry 
under local scalings of the probability density, derive the transformation law for the transition rates 
and interpret the thermodynamic force as a gauge potential. A widely accepted expression for 
the total entropy production of a system arises as the simplest gauge-invariant completion of the 
time derivative of Gibbs's entropy. We show that transition rates can be given a simple physical 
characterization in terms of locally-detailed-balanced heat reservoirs. It follows that Clausius's 
measure of irreversibility along a cyclic transformation is a geometric phase. In this picture, the 
gauge symmetry arises as the arbitrariness in the choice of a prior probability. Thermostatics 
depends on the information that is disposable to an observer; thermodynamics does not. 



Introduction. — Open systems subject to dissipation 
are usually modelled through markovian dynamics and 
further characterized by nonequilibrium thermodynamics 
PP. The link between dynamics and thermodynamics is 
the concept of thermodynamic force. In this work we as- 
sume that dynamics enjoys a gauge symmetry and show 
that the thermodynamical sector of the theory arises very 
naturally from the requirement that physical observables, 
in particular the entropy production, are gauge invariant. 
The driving force plays the role of gauge potential. 

A gauge theory has an internal symmetry whose action 
leaves all physical observables invariant [5] . Strictly speak- 
ing, any symmetry leaves observables invariant, there 
comprising the choice of units and reference frames [3]. 
However, conventionally one refers to gauge symmetries 
as to supplementary redundancies of the variables, which 
are not related to space-time symmetries: hence "inter- 
nal". One copy of the internal space where the gauge 
symmetry acts, is thought to be attached to each point i 
of space-time, or more generally of a base manifold. The 
peculiarity of gauge symmetries is that they act locally, 
that is, pointwise differently. The tool to compare nearby 
transformations is the connection, or gauge potential, A. 
Its transformation properties arc employed to adjust ob- 
servables {e.g. the Lagrangian in field theories) in order 
to guarantee the internal consistency of the theory. 



Crucial gauge-invariant observables can be built by cir- 
culating the connection along closed paths 7 over the man- 
ifold, obtaining the Wilson loops of the theory 



Win) = ^ 



V exp (J> A 



(1) 
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The trace over internal group-representation indices and 
the path-ordering operator V will be inessential in our 
theory. Wilson loops are the finite version of the curva- 
ture of a connection, telling how different paths between 
two given points carry different transformations. They 
play an ever more prominent role in our understanding 
of many phenomena besides field theories, from adiabatic 
phases in quantum mechanics [1] and quantum computa- 
tion schemes [S] to a proposal of quantum gravity [5] . 

According to a celebrated theory by J. Schnakenberg PQ, 
the macroscopic external forces (affinities) that maintain 
a nonequilibrium system into a dissipative steady state 
are defined as circuitations of the thermodynamic force. 
Schnakenberg claimed that affinities are the fundamen- 
tal observables which characterize nonequilibrium steady 
states, bringing as evidence that they satisfy Onsager's 
reciprocity relations in the regime where forces and cur- 
rents are linearly related. Andrieux and Gaspard [7] com- 
plemented his insight with a Green-Kubo relation for the 
linear response coefficients; the author proved that affini- 
ties serve as constraints for a formulation of the minimum 
entropy production principle near equilibrium [5]. 
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Schnakenberg's theory applies to markovian evolution 
of a normalized probability density over a discrete phase 
space, obeying the master equation 



j j 



iPji^-WjiPifj)] (2) 



where Wij is a positive transition rate from vertex j to 



vertex i, and Ji 



Vjipi is the probability cur- 



rent. In spite of space-time, the base manifold is the dis- 
crete phase space V consisting of a number \V\ of vertices. 
We define for later convenience an order relation i < j. 
Microstates are pairwise connected by a number \E\ of 
edges ij e E, which are assigned a conventional orienta- 
tion ij = i ■*- j. We assume that the graph G = (V,E) is 
connected and that forward and backward rates are non- 
null along all edges. The prominent role played by circui- 
tations in Schnakenberg's theory is an indicator that it is 
a gauge theory, with the driving force 



Aij = log^-. 



(3) 



playing the role of the gauge potential. To our knowledge, 
until very recently the geometric nature of the thermody- 
namic force was confined to the mathematical literature 
[9] and to work by Graham jlOj , with no explicit reference 
to gauge invariance. In the year of writing, for continu- 
ous diffusive processes the interpretation of the force as a 
gauge potential has been put forward by Feng and Wang 
[TTj . Sagawa and Hayakawa [12] made a proposal for a 
gauge potential connecting noncquilibrium steady states 
along slowly driven protocols; differently from Feng and 
Wang's, their connection has null curvature. They also 
observe that "the gauge symmetry does not seem to play 
any important role" . We fill the gap, taking an orthogonal 
approach: we do not assume the connection to be given, 
and derive it as the most natural candidate which guaran- 
tees covariance of the master equation; the appearence of a 
gauge potential is a byproduct of the symmetry, a conclu- 
sion which in a way parallels Abe and Kaneko's analysis 
of driven quantum equilibrium states |13j . 

Priors. — Prior to the derivation, we provide a phys- 
ical and slightly phylosophical motivation for the gauge 
symmetry, with a detour into continuous variables. In the 
"informationist" approach to statistical mechanics, whose 
forefather is Jaynes [M], the Gibbs-Shannon entropy 



S (P) = ~T,Pi l0 §^ 



(4) 



is a measure of the ignorance that an observer has about 
the state of the system. Maximization of the entropy, sub- 
ject to constraints according to whatever pieces of informa- 
tion the observer gains from measurement, produces the 
most plausible distribution given that the vertices of the 
system are a priori equally likely, according to Laplace's 
principle of insufficient reason. There is a source of subjec- 
tivity, which Jaynes accepted as physical, related to the 



choice of the observables one sets up to measure. How- 
ever, there is a second one which is in-built and which 
made Jaynes uneasy in his earlier writings. From |14j : 

Laplace's "Principle of Insufficient Reason" was an 
attempt to supply a criterion of choice [. . . ] How- 
ever, except in cases where there is an evident el- 
ement of symmetry that clearly renders the events 
"equally possible" , this assumption may appear just 
as arbitrary as any other that might be made. Fur- 
thermore, it has been very fertile in generating para- 
doxes in the case of continuously variable random 
quantities, since intuitive notions of "equally possi- 
ble" are altered by a change of variables. 

He then advised to replace Laplace's with the maximum 
entropy principle. So doing, he swept the dirt under the 
carpet, as the Shannon-Khintchin's set of axioms for the 
entropy include equiprobability 15 , and Shannon's mono- 
tonicity axiom [T3J p. 630] makes reference to it. 

Moreover, alleged paradoxes are found in the continuous 
variables case, where the differential entropy [TOJ Ch.9] 



s (p) = ~ dx p(x)\ogp(x) 

•J X 



(5) 



has been a source of dismay |17H19j , for it is not invariant 
under a change of variables. In fact, letting x' >->• x(x') be 
invertible with jacobian J(x') = det(dx/dx'), one has 



S(p) = S(p') + (log J), 



(6) 



where the probability measure and its density transform 
respectively as a volume form and as a volume density, 

p'(x')dx' = p(x)dx, p'(x') = J(x')p(x(x')) . (7) 

Related to this is the following riddle (from [20l Ch.8]): 
if we pick a number x between 1 and 10 at random, the 
probability that it is smaller than 5 is x /2; but if we pick 
x' at random between 1 and 100, the probability that it 
is smaller than 25 is How is it possible that picking 
either a number or its square aren't equally likely? The 
solution to this puzzle is to recognize that the choice of an 
arbitrary prior is congenital. It hides in that "at random" 
which is the continuous counterpart of Laplace's principle: 
in the first case we assume x is uniform, so that the prior 
is i/iodx; in the second we assume that x' = x 2 is uniform, 
with prior i/ioodx' = i/soxda;. Formally, in order to make 
eq.(j6]) mathematically sound, one will interpret p(x) as 
the Radon-Nicodym derivative of the probability measure 
p(x)dx with respect to the arbitrary prior da;. A change 
of variables corresponds to a change of prior. 

This is not, as Jaynes thought, an artifact of continu- 
ous variables. Think for example of a dice. Basing on 
visual impressions — which, by the way, are the result 
of a measure process — we might be tempted to assign 
equal probabilities !/6 to each face. However, if we knew 
that an incredibly huge mass was hidden near one of its 
corners, due to friction with air and the inelastic impact 
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with the gaming table, we would have sufficient reason 
to believe that the three faces which are adjacent to the 
loaded corner will have probability approximately l /3, and 
the others near zero. Our gambling strategy will depend 
on this prior knowledge. As a way out, in spite of invok- 
ing measure theory, up to additive constants we might just 
regard (neg)entropy as a special case of relative entropy 

s{p\\^ = Y J p^g Pl iP ( r\ (8) 

i 

with respect to a uniform prior p,- pr ^ = l^l -1 . We refer the 
reader to Banavar and Maritan's work |21j for some nice 
physical implications of working with relative entropy. 

The physical rationale is that the quantification of the 
entropy of a system depends on the choice of the under- 
lying degrees of freedom. If we assume that all configu- 
rations of positions and momenta of a number of classical 
particles are equiprobable, we implicitly coarse-grain the 
atomic and subatomic structure. The question "how much 
entropy is within a body" makes no sense on its own, since 
we can always go deeper into the inner structure of matter, 
according to the resolution of our "gedanken-apparatus" . 
Thermostatics depends on the prior. However — and here 
comes the key point — if we put a gas in contact with 
heat reservoirs, the process will occur in exactly the same 
manner, irrcgardless of our quantification of the system's 
entropy. Thermodynamics is independent of the prior. 

This is the gauge principle we assume. A gauge trans- 
formation is a change of priors. Gauge fixing means to 
choose a prior; it is analogous to the choice of a posi- 
tion with respect to which we measure displacements, with 
the important difference that in Newtonian mechanics the 
choice of a reference frame is absolute, in gauge theories 
the choice of a "reference frame" varies point by point. 

Gauge transformations. — With the r.h.s. of eq.© 
in mind, we postulate that 

Pi = e-"(H (9) 

is a symmetry of the theory, with ipi e R. For sake of 
consistency, we need to prescribe transformation laws for 
all of the objects which partecipate to the master equation. 

First, consider transition rates. We assume an edgewise 
and linear transformation law w[a = VijWije ipj , where we 
singled out e lpj , without loss of generality. Notice that the 
special case with all i?y = 1, namely 

w' l3 = w lJ e ipj (10) 

leaves the currents invariant, 

Jij = Jij, (11) 

The following graph-theoretical analysis proves that 
eq. (|10j) is the most general edgewise linear transformation 
law compatible with eq.©. 

Under the assumptions of connectedness and of nonva- 
nishing rates stated above, there exists a unique steady 



state p* of the master equation, which makes the r.h.s. 
of eq.© vanish. An explicit expression for p* , due to 
Kirkhhoff , is known [TU22] : we shall use it to constrain the 
Vij's. The recipe goes as follows, (i) Consider a rooted 
spanning tree Ti , i.e. a maximal set of edges with no cy- 
cles and with one preferred vertex i chosen as the root; all 
edges are oriented so that there exists exactly one directed 
path from any other vertex to i. (ii) Take the product of 
the transition rates along its edges, 

K w {T.i) = Yl w jk- 

jkeT, 

(iii) Form the sum over all possible rooted spanning trees 
to obtain the polynomial Zi = Y.T ^wiTi), which is ho- 
mogeneous of degree \V\ - 1 in the transition rates, (iv) 
The steady state is p* = Zi/ Y,j Zj. By construction, each 
vertex of the graph, but i, is the starting point of exactly 
one edge of Tj , so that we can factorize 

Z , i = e*-**'£ir v (T i )ir v ,(Td, (12) 

Ti 

where 4> = T,jtv Vj ■ By eq. © , Z[ must be proportional to 
the analogous expression for e~ Vi Zi . We obtain 

EMTO = c 1 -^Y,^v(T i )n w (Ti) (13) 

Ti Ti 

where c is a proportionality constant. Since all transition 
rates are positive and the equality must hold Vw's , if fol- 
lows that ir v (Ti) = c^^ 1 independently of the spanning 
tree. Furthermore, the transformation law should be uni- 
versal, i.e. it should not depend on the specific graph. 
As graphs become larger, tipically the number of span- 
ning trees grows exponentially in the number of edges of 
the graph — whereby "tipically" loosely means "for most 
graphs" 23J. This entails that the number of equations 
specified by -K v (Ti) = c'^' -1 becomes enormously larger 
than the number of the unknowns. The only universal 
solution is My = c, Vr/; the constant can then be scaled to 
unity with a redefinition of the time unit r -*■ t/c. □ 

We now face a seeming paradox. In fact, considering the 
transformed master equation p[ = Y,j J[j an d keeping into 
account eq.® and eq. (|llj) . we obtain an equation which 
is not equivalent to the starting master equation, eq.©! 
The solution delves into the geometrical interpretation of 
summation symbols. We introduce the incidence matrix 

(+1, if j < k, k = i 
-1, i£j<k,j = i , (14) 
0, elsewhere 

and rewrite the master equation as a continuity equation 

p+ x dJ = 0. (15) 

Technically speaking, the incidence matrix is a boundary 
operator which maps edges into their boundary vertices. 
Normalization of the probability can be written as 

5>=°0/>=l (16) 

ieV 
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where we introduced one further boundary operator d = 
(1,1,. ..,1), which maps vertices to the connected com- 
ponent of the graph they belong to. Although this latter 
definition is rather trite, it gets more interesting when the 
graph has several disconnected components. Notice that 
°d l d = 0, from which conservation of probability follows. 
Therefore "summation over i" has different geometrical 
meanings according to the context. 

Strictly speaking, pi should not be considered as a num- 
ber, but rather as a one-component vector which lives in 
the internal vector space = R which is attached to ver- 
tex i. The gauge transformation cq.® is interpreted as 
a linear change of basis in It follows that we should 
consider the boundary operator's entries as linear maps 
on ^j, which therefore transform according to 

% = e~ Vi % % = %. (17) 

With this prescription, eq.s (TT5|) and (|16l) are covariant. 
To simplify the notation a bit, we introduce a modified 
sum symbol such that 

E'=E^ (is) 

i i 

This modified symbol is crucial for the up-coming result, 
so let us further linger on it. Consider the average of a 
gauge- invariant vertex function f = f (a scalar field), 

i i 

Requiring gauge invariance V/ yields the transformation 
law for the summation symbol. In other words, while the 
probability measure ( • ) p is gauge invariant, the proba- 
bility density pi is not, in analogy with the continuous 
variables case, see eq.©. 

We now focus on the Gibbs-Shannon entropy, which 
transforms according to 

SS = S'[p'(r)] - S[p(r)] = (^> p = -S(p' \\ p) (20) 

where S' is calculated using Y! ■ On the r.h.s., the trans- 
formation law is succinctly expressed in terms of relative 
entropy. Remarkably, while relative entropy is not a differ- 
ence of entropies, in this context it is naturally interpreted 
as (minus) the entropy change after a gauge transforma- 
tion. The rate at which the entropy of the system changes 
is subject to 

6& = ZM<Pi-<Pi)- (21) 

i<j 

In gauge theories, non-gauge invariant terms are adjusted 
with the introduction of a connection, which is an anti- 
symmetric edge variable Ajj = —Aji such that 

SAij = tpj ~(pi. (22) 

Once a connection is given, the tern0 

* = E-M« ( 23 ) 

i<j 

x In this case, the sum has meaning of a bilinear form from the 
space of edges to real numbers, so no gauge transformation is needed. 



has a transformation law which balances eq. ([21[) . mak- 
ing S + a invariant. In principle, connections can be con- 
structed as convex linear combinations of terms such as 

log 4, log^, (24) 
P* w< 

where u>i = Y,k w ki is the average frequency of a jump out 
of site i. So, for example, adding Y,%<j Jij\°gPi/Pj yields 
the relative entropy with respect to the steady state. The 
latter plays an important role in the theory of Markov 
processes as a Lyapunov functional 1 Sec. V]; fitly, it is 
gauge invariant, while entropy per se is not. However, the 
options listed above are, technically speaking, exact: they 
are differences of vertex functions, so that their circui- 
tations vanish, thus making the graph's geometry rather 
dull. As a further consequence, gauge invariant terms ob- 
tained this way vanish at the steady state. 

A good candidate as a "truly edge" connection variable 
is given by the driving force, defined in eq.©- Although it 
is not the only antisymmetric edge variable that one could 
engineer which transforms according to eq. (|22|) . it is cer- 
tainly the simplest. Then S + a coincides with Schnaken- 
berg's total entropy production [TJ eq.(7.6)], 

otot = S + a = X-A,luo-^i. (25) 

which is widely accepted as the entropy production rate 
of a Markov process [241125] . In this setting a arises as the 
simplest term which completes S into a gauge invariant 
quantity and which does not vanish at the steady state. 

A gauge transformation will result in a shift of a total 
time derivative from a to S, with a consequent redefinition 
of the internal entropy and of the entropy flow towards the 
environment. For example, letting tpi = \ogp*, we obtain 

S' = -S(p\\p*), a'=EJylog^4, (26) 

whose microscopic analogues along single stochastic tra- 
jectories have been interpreted by Esposito and Van den 
Broeck as non-adiabatic and adiabatic terms, obeying de- 
tailed fluctuation theorems |26j . Gauge transformations 
of fluctuation theorems will be discussed in a later work. 

Wilson loops. — From a geometrical viewpoint [27], 
not only A provides a connection over the manifold, but it 
also constitutes a measure of the oriented lenght of paths 
along chains of edges rj = (i n i n -i, ■ ■ ■ ,ii«o) 5 

Efo) = t A ^ = f A - (27) 

Since the lenght is additive upon composition of paths, 
the real positive numbers obtained by exponentiating £(??) 
can be thought of as elements in the multiplicative group 
of real positive numbers (R + , x), which is the gauge group 
of the theory. In the representation theory of groups, 
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group elements are not seen as "static" objects, but rather 
as "active" linear maps; they act on vectors ipi which live 
in the internal vector spaces \&$. Such vectors acquire 
phases as they are parallel transported along paths, thus 
connecting far-apart vertices, 

V>i„ = exp£(r?) ip io , (28) 

where tpi n represents the result of parallel transport along 
path r\. In our case, due to the very simple gauge group, 
the displaced vector is just a real number and parallel 
transport produces a scaling factor. The intepretation of 
group elements as linear maps further entails that new 
equivalent representations can be obtained by perform- 
ing basis transformations in one per each vertex: this 
yields a gauge transformation. In the case at hand, such a 
basis change amounts to an orientation-preserving rescal- 
ing %p[ = e~ lf>i Tpi. Transformed vectors are parallel trans- 
ported according to ip'i n = exp ip' ig , where £' is a new 
representation of the group element, defined in terms of 
a transformed connection A'. Requiring equivalence with 
eq.(|28|) for any possible path rj yields the transformation 
law for the vector potential, eq. (|2"2"|) . Grossly, this intro- 
duces the geometrical framework for gauge theories. 

Gauge transformations define an equivalence relation 
" ~ " between gauge potentials; so, for example, the adi- 
abatic force log(wij p* /utjip*) [26] is gauge-equivalent to 
ours, eq.Q. The connection is said to be exact when it 
is equivalent to A'^ =0. It is well known [5 HTT1|2"8] that 
equilibrium systems are characterized by an exact poten- 
tial. In fact, when Ay = ipi -<fij, the steady solution of the 
master equation is p* ex e~ Vi , as direct substitution into 
eq.® shows. Detailed balance follows: 



Along closed cycles 7, with io = i n , the exponentiated 
lenght is a Wilson loop. When Wilson loops are all unity 
the connection is exact and the oriented lenght of an open 
path only depends on the extremal vertices, and not on 
the particular path which connects them, for which rea- 
son the connection is said to be flat. As a remarkable 
consequence, Kolmogorov's criterion |22[|2"§] is equivalent 
to all Wilson loops being equal to unity. Hence detailed 
balanced systems can be seen as the special class of mod- 
els with a flat connection, with zero curvature; they all 
belong to the same equivalence class. 

On a discrete state space, knowledge of a finite num- 
ber of Wilson loops suffices to characterize the connection. 
The so-called Mandelstam identity 

W{ 11 o l2 ) = W{ ll )W{ l2 ) (30) 

allows to compose loops. A basis of loops can be found 
this way. Consider an arbitrary spanning tree T of the 
graph — this time with no preferred root and orientation. 
Let i a j a be one of the edges which do not belong to T. 



By definition, adding i a j a to the spanning tree generates 
a cycle j a , which can be oriented according to the orien- 
tation of i a j a - By Euler's formula, there are \E\ - |V| + 1 
such cycles. It is a basic graph-theoretical result that any 
loop can be decomposed in terms of the 7 Q 's [3TJ Part I, 
Ch.4]. Let efj be +1 if ij = i a j a , _ 1 if ji = iaja, otherwise 
it is zero. It can be shown that 

log YlW{ la )^ ~ A^. (31) 

a 

Hence Wilson loops allow to reconstruct the gauge poten- 
tial, up to gauge transformations 30 . By eq. (l3T1) . the 
choice of a spanning tree fixes the gauge by selecting one 
particular representative in the equivalence class of Ay . 

Spanning trees also allow to give a physical interpreta- 
tion of the connection, as follows. Any graph which co- 
incides with a spanning tree, E = T, has no cycles, hence 
it can only accomodate equilibrium systems. Then there 
exists a vertex function ipi = fiiii such that 

Wij/wji = e^-«*>, ijcT, (32) 

where we introduced an inverse temperature j3, in units 
of Boltzmann's constant. The inverse temperature and 
the energy Ui are determined up to an energy shift and a 
rescaling of units, Ui -* k{ui + «), /3 -*■ AT 1 /?. In general, 
adding further edges i a ja to the graph will not result in a 
detailed balanced system, unless we fine-tune their rates. 
We then define a new set of temperatures (3 a , such that 

m a jjw jaia = e Mu ^- u ^\ iaja e E\T. (33) 

We just proved that the thermodynamics of any collection 
of transition rates can be described in terms of at most 
\E\ - \V\ reservoirs, each at its own temperature, satisfy- 
ing the condition of local detailed balance [32 . In this 
"minimal" case each transition is due to the interaction 
with exactly one reservoir. This ansatz allows to recast 
the basis Wilson loops in this form 

W{ la ) = exp[(/3-/3 Q )(^-^J]. (34) 

Therefore, temperature differences are the fundamental 
thermodynamic forces of nonequilibrium systems, as one 
could expect. Since there is no external time-dependent 
driving, which would result in time-dependent transition 
rates, no work is performed by an external agent along 
one single realization of the process, and by the first law 
of stochastic thermodynamics [24], along a transition the 
energy gap Su coincides with the heat exchanged 5q. It is 
then illuminating to rewrite the geometric phase as 

IogW( 7 ) = jf^, (35) 

yielding Clausius's measure of irreversibility along one re- 
alization of a cyclic irreversible process. The lenght E(^) 
is the entropy exchanged with the environment along any 
trajectory which performs a sequence of jumps, whichever 
the jumping times might be. This notion is completely in- 
dependent of the time parametrization of the trajectory: 
it is purely geometrical. 



p-5 



Matteo Polettini E-mail: polettiniObo . infn. it 



Final remarks. — The above construction can be eas- 
ily generalized to Markov processes with time-dependent 
transition rates and to time-dependent gauge transforma- 
tions. In this respect, our formalism has evident points 
of contact with stochastic pumping along cyclic protocols, 
whose geometrical nature has been recently studied 33 . 
It would be a conceptual advance to give a unified descrip- 
tion of both aspects of NESM. We notice in passing that 
Sinitsyn 34, §6] makes a remark on gauge transformations 
applied to the current generating function, arguing that 
they follow from the modification of the "prior" currents 
which have flown before a given initial time. 

Regarding the nature of gauge transformations, for con- 
tinuous variables they have been shown to follow from co- 
ordinate changes. Thus the gauge group could be seen as 
(a subgroup of) the group of diffeomorphisms. It has been 
a matter of disagreement [35] [6] §2.1.3] whether diffeomor- 
phisms and gauge transformations should be considered by 
the same standards; the diatribe mainly revolves around 
gravity and its formulation as a local affine theory [36] . 
The identification of ours as a gauge transformations is 
justified by the usage of the gauge machinery, which is 
analogous to well-established practice for the formaliza- 
tion of geometric phases in QM and of electromagnetism 
as a U(l) gauge theory. Employing analogies with the 
latter, C. Timm [37] discussed a slightly different gauge- 
theoretic structure for master equations. 

To conclude, while we are conscious that the very sim- 
ple gauge group makes the geometrization of irreversible 
thermodynamics unnecessary for all practical purposes, it 
allows to better appreciate the importance of macroscopic 
affinities as fundamental observables, and it might serve 
as a good starting ground for later generalizations. We 
point out that a Schnakenberg-type analysis is still lack- 
ing for quantum nonequilibrium systems, either described 
by a Lindblad-type equation or a by a more general in- 
teraction of a system with reservoirs of quantum degrees 
of freedom. It is tenable that excursions to the quantum 
world might require more interesting gauge groups and a 
more pertinent application of gauge theory. 



Joint work with Massimiliano Esposito is at the basis 
of the physical interpretation of loops; the author is also 
grateful for his useful comments. Stimulating discussion 
with A. Vulpiani and A. Montakhab is recognized. 
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